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We update and compare the capabilities of the purely leptonic mode B ^ tv and the semileptonic mode B — > Dtv 
in the search for a charged Higgs boson. 



- .- 1. INTRODUCTION 

oo ; 
o . 

■ Supersymmetric extensions of the standard model (SM) - or more generally extensions that require the existence 
, of at least one additional Higgs doublet - generate new flavour-changing interactions already at tree-level via the 
^ ' exchange of a charged Higgs boson. The coupling of to fermions grows with the fermion mass. It is thus natural 
to look at (semi)leptonic B decays with a t in the final state to try to uncover this type of effects. In a two-Higgs- 
I doublet-model (2HDM) of type II, where up-type quarks get their mass from one of the two Higgs doublets and 
fS| ' down-type quarks from the other one, effects are entirely parametrized by the H'^ mass, Mh, and the ratio 
of the two Higgs vacuum expectation values, tan/3 = Vu/vd- They can compete with the exchange of a W'^ boson 
<~i for large values of tan/3 [1]. In the minimal supersymmetric extension of the SM (MSSM), the tree- level type-II 
structure is spoilt by radiative corrections involving supersymmetry-breaking terms. The effective scalar coupling gs 
Oh then exhibits an additional dependence on sparticle mass parameters when tan/3 is large [q = w, c) [2, 3]: 



ttH+ o /7^r< T/ '^brrir r_ , i r_ i , , M| tan^ /3 1 

iJ,,, = -2V2G,V,,^gs [Uhn] [ruM + h.c, gs = (1 -f So tan/3)(l + tan/3) ' 

[ where eo,r denote sparticle loop factors. The correction induced can be of order one. However, the access to the 
' Higgs sector remains exceptionally clean. In Eq.(l), gs has been normalized such that it gives the fraction of effects 
in the B ^ TV amplitude, which is very sensitive to exchange: B{B — > tv)/B{B tv)^'^^ = |1 — gs]"^. The 
B Dtv channel is less sensitive (though better in this respect than other modes such as i? ^ D*tv) but, as we will 
\ see, exhibits a number of features that make it, too, play an important part in the hunt for the charged Higgs boson. 

o ■ 

> . 2. B{B Dtu) versus B{B tu) 

' The current capabilities of B{B Dtv) and B{B tv) to constrain effects are compared in Fig.l for 
5S > (as is typically the case in the MSSM or the 2HDM-II). The lower sensitivity of the B — > Dtv mode 
comes from the different momentum dependence of the Higgs contribution with respect to the longitudinal one^ : 
{dT{B Dtv) / dq'^)^ii '^^ cx |1 — gs{q^ /Mg)/{1 — mc/m},)\'^ with q = pB — Pd- On the other hand, the theory 
prediction for B{B —> tv) suffers from large parametric uncertainties from the CKM matrix element Vut and the 
B decay constant /b- In contrast, Vet is known with better than 2% accuracy from inclusive B Xc£v {£ = e,^) 
decays, \Vcb\ = (41.6 ± 0.6) x 10~^ [4], and the form factors f+{q^) and fo{q^) describing the B ^ D transition are 
very well under control, as we now discuss in more detail. 

To this end, we introduce the following conformal transformation: 



, 2 ViMB + Mp)^ - g2 _ ^(Mb + Mp)^ - to 

z(q ,to) = — , = , (zl 

^{Mb + Mp)^ -q^ + ^{Mb + Mp)^ - to 



^Note that the latter, though holicity-suppressed, is still (slightly) larger than the transverse contribution for all values. 
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Figure 1: Left: B{B tu) as a function of gs- Light gray band: 8"""^ = (1.51 ±0.33) x 10"'' [5]. Gray (blue) band: S"' for 
\Vub\ ~ (3.95 ±0.35) X 10~^ [4] and fs = 216±38MeV [6] (in this last case, we add errors linearly to stay on the conservative 
side). Right: R = B{B -> Dtu)/B{B D£u) as a function of gs- Light gray band: = (41.6 ± 11.7 ± 5.2)% [7]. Dark 

gray (dark blue) band: i?"' with the HFAG vector form factor before 1CHEP08 [8]. Gray (blue) band: T?**" with the HFAG 
vector form factor after ICHEP08 [5]. The dashed lines indicate the 2 and 3-sigma limits. The ratio mc/mt in the MS scheme 
has recently been determined with very high accuracy: mc/mt — 0.2211 ± 0.0044 [9]. We inflated the error on this number 
and set mc/mi, = 0.22 ± 0.01 to reduce the discrepancy with the HFAG estimation [5]. 



which maps the complex plane, cut along > {Mb + Md)^ , onto the disk \z\ < 1. The form factors /+ and /o 
are analytic in z in this domain, up to a few subthreshold poles, and can thus be written as a power series in z after 
these poles are factored out {i — +,0) [10]: 

M,^) - Hito) + alito)ziqMo) + ...] , (3) 

where the function Pi gathers the pole singularities and an arbitrary analytic function can be factored out as well. 

This parametrization has been used in Ref.[ll] with the choice to = g^^x = {Mb — Md)'^ , together with heavy-quark 
spin symmetry inputs, to derive the following ansatz for the vector form factor: 

U{q') ^ M^±MRv^{q'), V,{q^) = 5(1) [l - Sp'ziq^to) + (Slp^ _ lQ)z[q\ tof - (252^^ _ 84)z(g2, tof] , (4) 

where Vi is defined such that it reduces to the Isgur-Wise function in the heavy-quark limit and Q{1) = Vi{q^^^). The 
parameters |Vcfc|5(l) and can be determined from B Dlv experimental data. Before this summer, the HFAG 
averages [8] based on BELLE, CLEO, and ALEPH data read: \Vab\G{l) = (42.3 ± 4.5) x lO'^ and = 1.17 ± 0.18 
(with a \Vcb\Qil)-p'^ correlation of 0.93). The recent BABAR resuhs [12] and [13] have now been included, leading 
to a substantial improvement [5]: \Vcb\G{l) = (42.4 ± 0.7 ± 1.4) x 10"^ and ^ 1.19±0.04±0.04 (with \Vcb\G (l)- p'^ 
correlation 0.88). The old and new vector form factors are compared in Fig. 2 (left), where we have defined as usual 
w; = (M| + Ml, - q^)/(2MBMD). 

For the scalar form factor, we adopt the ansatz of Ref. [14]: 

Uq^) ^ i!fi±il^^^5i(g^) = , , , , , ^ . K(to) +a?(io).(g^to)] , (5) 

Mb + Md z{q^,M{)z{q^,M^)(po(q^,to) ^ 



where to = {Mb+MdY [l - \J\ - {Mb - MdY/{Mb + MdY j such that [^Imax is minimized, Mi = 6.700 GeV and 
M2 = 7.108 GeV [15] are the subthreshold poles, and 4>o is obtained from Eq.(lO) of Ref.[14] setting = and rj = 2: 

, ^^_J 2{MI-Myf ^{MB+Mof-q^ z{q\Qf ( z{q\to) \-''^ ( zjq^jMB - MnY) \-'" 
^o[q,to)-y {{MB + Mnr-toV/^ {q^Y \to-q') \{MB^MnY^q^ ) ' 

Following [16], we truncate the series (3) after the first two terms. This is motivated by the fact that |z|max = 
0.032 and that a similar parametrization for /+, when fitted to experimental data, produces the same result as 
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Figure 2: Vector (left) and scalar (right) form factors corresponding to the |Vcb|t/(l) and determinations of HFAG before 
(dark gray/dark blue)[8] and after (gray/blue) [5] ICHEP08. With the new determination, the errors on |Vc6|Vi and IVctlS*! 
are smaller than 4% and 7%, respectively. 




Figure 3: 95% C.L. exclusion zones in the (MH,tan/3) plane from B{B tv) (dark gray/dark blue) and R (gray/blue) in a 
2HDM (left) and in the MSSM with eo — 0.01 and e,- — (right). The exclusion limits are directly read from the gray (blue) 
bands in Fig.l. They differ from those usually found in the literature in that experimental and theory errors are not simply 
added in quadrature and the dependence of the errors on the contribution is taken into account. 



Eq.(4) in very good approximation [16]. Then, |Vc6|aQ(to) and |Vc6|a5(to) are determined imposing the conditions (i) 
\Vcb\Si{0) = \Vcb\Vi{0) and (ii) IKbl^iteax) = (4.24 ± 0.27)% (corresponding to l^hl = (41.6 ± 0.6) x 10-^ from 
B Xciv [4] and Si^q^^^ = 1.02 ± 0.05 from HQET [16]). The scalar form factors obtained in this way from the 
old and new |Vcb|Vi are not very different, as one can see on Fig. 2 (right). 

The recent progress on |Vct|X^i allows to reduce the errors on the SM predictions for the two B Dtv branching 
fractions: B{B- D^t-v)^^'^ = (0-70t°:^^) %, ^(5° ^ D+r-vf^^ = (0.65t°:g^) % (differing essentially due to 
tbo 7^ Tg+). The errors from [KblX^i, however, already cancel to a large extent in the ratio R = B{B — > Dtv)/B{B 
Dliy), which is why the nice improvement in Fig. 2 has little impact on Fig.l (right), already dominated by the error 
on S'i(g'max)- R^^^ — 0.31 ± 0.02. This estimation is compatible with the one obtained from lattice methods: 
Rfjf^ = 0.28 ± 0.02 [17]. Note that replacing condition (ii) by a constraint on Si{ql^^J/Vi{ql^^J from HQET would 
lead to a similar error on R. Still, an interesting 95% C.L. bound on gs can already be obtained from R: gs < 1.79, 
complementary to the bounds from B{B — > ti/): gs < 0.36 U 1.64 < gs < 2.73. The corresponding exclusion zones 
in the (M^f,tan/3) plane are depicted in Fig. 3. The error assigned to Si{q^^^) is quite conservative, so the above 
constraints are robust. At the three-sigma level, it is not possible to extract any interesting bound from R yet, but 
its experimental knowledge is expected to improve in the near future. Its role to constrain effects will then of 
course depend on the new central value. For the moment, a 15% measurement with the same central value would 
exclude gs > 0.29 at the 95% C.L.. 
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Figure 4: dr{B ~* DTu)/dw (left) and dr{B Dvt[-^ -kv]) / dETrdcosOD-ndw with = 1.8 GeV and cosOdt, = -1 (right) 
for gs = Q (gray/red), gs = 0.35 (light gray/light blue), and gs ~ 1.75 (dark gray/dark blue). These values are still allowed 
by B{B — > Dtv) and B{B tv) at the 95% C.L.. The various curves have been obtained using the more recent HFAG vector 
form factor [5]. The lighter bands take all errors into account, while the darker bands only take into account the error on 
Si ((/max)- One could of course also normalize the above differential distributions to dV{B — > D(.v)/dw to reduce the impact of 
the errors on /+. 

3. B ^ Dti^ differential DISTRIBUTIONS 

If a hint for a charged Higgs boson is seen at the branching fraction level, B Dtv has a great advantage over 
B —f TV. it allows to analyze the same data points on a differential basis, better suited to discriminate between 
effective scalar-type interactions and other effects. The dT{B DTv)/dq^ distribution, in particular, has already 
been studied in great detail [18]. The polarization of the t is also known as a analyzer [19], yet it requires the 
knowledge of the r momentum, which cannot be accessed at B factories as the r does not travel far enough for a 
displaced vertex and decays into at least one more neutrino. 

A straightforward way to nevertheless exploit the sensitivity of the r polarization to effects and at the same 
time retain the information from the spectrum is to look at the subsequent decay of the r into a pion and a 
neutrino [16]. The direction of the pion is indeed directly correlated with the polarization of the t. Integrating over 
the neutrino momenta, we end up with a triple differential decay distribution dT{B — > Dvt[—^ ttu]) / dq'^ dE^^d cos 0dtt ■ 
An explicit formula is given in Eqs. (9-11) of Ref. [16] (with = /+ audits = /o). Its sensitivity to 55 is illustrated in 
Fig. 4 for i?7r = 1.8 GeV and cos9]Jt:^ = — 1. For comparison, we also display the q^ spectra corresponding to the same 
gs values. Of course, in practice, one should not fix i?^ or Odtt, but rather perform a (unbinned) maximum likelihood 
fit of the triple differential decay distribution to the available data points. The information from the q^ spectrum in 
the dominant r Ivv decay channel should also be included in the fit to make the most out of experimental data. 

4. CONCLUSION 

The form factors f+{q^) and fo{q'^) in the B — > Dtv transition are under good control. As a result, the ratio 
R = B{B -> Dtv)/B{B -> Div) can be predicted with 7% accuracy in the SM: R^"' = 0.31 ± 0.02, where the 
5% uncertainty on the scalar form factor at zero recoil >S'i((?max) is the main error source. This allows to derive 
useful constraints on the effective coupling gs- Together with the constraints from B{B tv), we obtain: 
gs < 0.36 U 1.64 < gs < 1.79, i.e., the window around gs = 2 left over by B{B tv) is now nearly completely 
excluded by R alone. These bounds should be strengthened soon thanks to the current considerable experimental 
efforts on both modes. In this respect, one should pay particular attention to the B —>■ Dtv differential distributions 
as these are especially well-suited to discriminate between effective scalar interactions and other types of effects and, 
if the former are seen, to extract the coupling gs with good precision. 
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